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Abstract 

C^J The recent discovery of free-floating planets and their theoretical in- 

terpretation as celestial bodies, either condensed independently or ejected 
from parent stars in tight clusters, introduced an intriguing possibility. 
Namely, that exoplanets are not condensed from the protoplanetary disk 
of their parent star. In this novel scenario a free-floating planet inter- 
acts with an already existing planetary system, created in a tight cluster, 
and is captured as a new planet. In the present work we study this in- 
teraction process by integrating trajectories of planet-sized bodies, which 
encounter a binary system consisting of a Jupiter-sized planet revolving 
around a Sun-like star. To simplify the problem we assume coplanar orbits 
for the bound and the free-floating planet and an initially parabolic orbit 
for the free-floating planet. By calculating the uncertainty exponent, a 
quantity that measures the dependence of the final state of the system 
on small changes of the initial conditions, we show that the interaction 
process is a fractal classical scattering. The uncertainty exponent is in 
the range (0.2 — 0.3) and is a decreasing function of time. In this way we 
see that the statistical approach we follow to tackle the problem is justi- 
fied. The possible final outcomes of this interaction are only four, namely 
flyby, planet exchange, capture or disruption. We give the probability of 
each outcome as a function of the incoming planet's mass. We find that 
the probability of exchange or capture (in prograde as well as retrograde 
orbits and for very long times) is non-negligible, a fact that might explain 
the possible future observations of planetary systems with orbits that are 
either retrograde or tight and highly eccentric. 

Keywords: planetary systems, numerical methods, statistical 
methods 
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1 Introduction 



The three-body problem has been a topic of active research since the seventeenth 
century, initially most of the effort focusing on bound solutions and their long 
time stability. This problem was "solved" through the work of (a) Poincare 
1893, who showed that the general solutions of even the simplest variant of the 
problem (apart from highly symmetric configurations), the restricted planar 
circular three-body problem, cannot be written in the form of convergent series, 
and (b) Marchal and Bozis (1982), who showed that in the general three body 
problem zero velocity curves do not exist, so that one of the bodies can always 
escape to infinity. However the "inverse" problem, i.e. the interaction of a body 
coming from "infinity" with a binary pair, has attracted less attention. From 
the astrophysical point of view, there are three distinctive special forms of this 
problem: (a) the binary members and the incoming body have similar masses 
(interaction of a "passing" star with a binary star), (b) the binary consists of 
a massive and a small body and the incoming body is massive (a passing star 
interacting with a planetary system) and (c) the binary consists of a massive 
and a small body and the incoming body is small (a passing planet interacting 
with a planetary system). The first two forms have already been discussed in 
the literature (e.g. see Boyd & MacMillan, 1993; Mikkola, 1994; Donnison, 2008 
for the first, Astakhov and Farrelly, 2004 for the second). However the third one 
has attracted limited interest up to now (e.g. see Donnison, 2006), probably 
because planets were considered always to be bodies revolving about stars. The 
discovery of free floating planets (FFPs) a few years ago changed this picture. 

The initial report by Zapatcro Osorio ct al. (2000) on the observation of free- 
floating planets with masses mj < m < 3m j (where to,/ is Jupiter's mass) was 
regarded initially with scepticism, since it was not consistent with the generally 
accepted definition of a planet. However subsequent microlensing observations 
(Sumi et al., 2011) showed that celestial bodies with masses to m mj really 
do exist in the Milky Way, and for that in numbers exceeding those of its 
stars, making the interaction of FFPs with already existing planetary systems 
a noteworthy process. In any case, recent publications (e.g. see Adams et 
al., 2006; Malmberg it et al., 2007) have shown that the probability of such an 
interaction is non-trivial within clusters of newborn stars, where a planet ejected 
from a young planetary system (or condensed independently) may remain within 
the cluster for long time intervals. In the present paper we study numerically 
the simplest possible model of such an interaction, namely a bound planet (BP) 
on circular orbit around a star, scattering a FFP approaching from "infinity" 
on a coplanar orbit. 

The interaction of a FFP with an existing star-planet binary system has only 
four possible outcome states: the incoming body (a) goes to infinity, probably 
exciting the planet (flyby) , (b) replaces in orbit around the star the BP, which in 
turn goes to infinity (exchange) (c) is captured by the star-planet binary system, 
forming a two-planet planetary system (temporary capture) or (d) the system 
is disrupted (disruption) (Fig. [T]). If the initial total energy of the system is 
negative, as in our numerical experiments, the fourth outcome is ruled out and 
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Figure 1: Characteristic orbits of the scattering between a star - planet system 
and a free-floating planet (FFP) (FFP: light gray, bound planet (BP): dark 
gray, star: black). Flyby: a simple deflection of the FFP, exchange: the FFP 
displaces the BP and takes its place orbiting the star, temporary capture: the 
two bodies are temporarily bound to the star, disruption: the system breaks 
up and the three bodies move to infinity. Note that disruption is not possible 
when the system's center of mass energy is negative, which is the case in our 
numerical experiments. 



we are left with only three possible outcomes. 

From the purely dynamical point of view, however, the situation is more 
complicated. The interaction can be distinguished into "single-" and "multiple- 
encounter" events. In a single-encounter event the FFP interacts only briefly 
with the BP (e.g. for a time interval less than one unperturbed period of the 
BP) while in a multiple encounter event the FFP enters in orbit around the star 
and interacts for an extended period of time with the BP (typically much more 
than one unperturbed period). Therefore a single encounter event can result 
only in flyby or exchange, which may be thought of as "direct" outcomes. To the 
contrary, a multiple encounter event (temporary capture) is not a "definitive" 
outcome, since after some time interval one of the two planets may escape to 
infinity (e.g. see Marchal & Bozis, 1982). In other words, a temporary capture 
may ultimately decay to a flyby or an exchange, which, in this sense, may 
be considered as "indirect" outcomes. Therefore finally a more comprehensive 
classification of the possible outcomes may read as follows: flyby (direct and 
indirect), exchange (direct and indirect) and "long-time" capture. It is not 
possible to guarantee a permanent capture since, as we mentioned already, in 
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the general three-body problem closed zero-velocity curves do not exist (except 
from cases that are Hill stable). But, as we show in section 4, the statistical 
analysis of our results indicates that a percentage of temporary captures end 
up in "deep" captures, in the sense that the resulting two-planet system may 
remain stable for very long times. In fact, if a fourth body exists nearby or 
other forces (e.g. tides) are taken into account, the 3-body system might be 
stabilized. 

The paper is organized as follows. In the next section we describe the setting 
of our numerical experiments. In section 3 we show that the scattering process is 
a fractal one, we calculate the uncertainty exponent and we give the probability 
of each one of the possible outcomes of the scattering. In section 4 we present 
our main results, concerning (a) the dependence of the scattering outcomes on 
the mass of the incoming FFP, (b) the physics behind capture and (c) the orbital 
distribution (in the cases of capture and exchange) as a function of the mass 
and the impact parameter of the FFP. In section 5 we deal with the decay of 
temporary captures and in section 6 we find the initial conditions that lead to 
prograde or retrograde orbits. Finally in section 7 we discuss our results. 

2 Numerical experiments setup 
2.1 Initial configuration 

We perform numerical experiments of the three body scattering between a star- 
planct system and an incoming FFP. The three bodies are considered as point 
masses moving on coplanar orbits and are: a star of mass M, a BP of mass mj 
in orbit around the star, and the FFP with mass m similar to that of the BP, 
incoming from infinity. The initial configuration of the scattering experiment is 
shown in Fig. [2] The star is initially placed at rest at the origin of the inertial 
reference frame Oxy and the BP is placed at a distance ro with initial velocity, 
vq, corresponding to a circular orbit around the star. The FFP is coming from 
a fixed distance of 40ro with impact parameter d and initial parabolic velocity v 
with respect to the star. This distance is sufficient to ensure that the FFP and 
the original planetary system are not initially interacting strongly. The initial 
velocity, v, of the FFP is estimated as a function of d through the relationship: 



where fi = (mM)/ (m + M). Keeping all other parameters fixed, the initial state 
of the system is determined by the impact parameter of the FFP, d, and the 
initial phase of the BP, <fi. Note that d can take both positive and negative values. 
By convention, we assign a negative value to d when the angular momentums of 
the BP and the FFP have different signs (y > in Fig. [5]) and positive otherwise 
(y < in Fig. [2|. The initial values of the impact parameter are restricted in 
the interval — 7ro < d < 7ro, where exchange and temporary capture are mainly 
expected to occur (see e.g. Donnison, 1984a, b). Within the configuration of 
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Figure 2: Initial configuration of the scattering between a star-planet system 
and a FFP. The figure shows the initial position of the three bodies with respect 
to an inertial reference frame Oxy. The BP is initially on a circular orbit around 
the star and the FFP is coming with parabolic co-planar velocity from a fixed 
distance, sufficient to ensure that it is not initially interacting strongly with the 
planetary system. The initial state of the system is determined by the impact 
parameter, d, of the FFP and the initial phase, </>, of the BP. Note that d is 
negative for retrogade orbits (y > 0) and positive for prograde orbits (y < 0). 



Fig. 2 the total energy of the system at the center of mass is always negative, 
although it depends weakly on <fr and d. 

It should be noted that the (dimensionless) interval of impact parameters 
chosen here is the one that yields dynamically interesting results. FFPs passing 
at a larger distance have minor effects on the BP. That is why the most promising 
site for this kind of interaction is in dispersing clusters (see, e.g., Perets & 
Kouwcnhovcn, 2012). But the dimensional value of impact parameters in AUs 
depends on ro- If one is taking as example our solar system, then he could 
consider as "bound planet" (BP) anyone of the giant planets. When considering 
Jupiter as BP the impact parameter is 35 AU, while when considering Neptune 
it is 210 AU. 



2.2 Equations of motion and integration 

Our dynamical system is an autonomous one with six degrees of freedom and is 
described by the Hamiltonian: 



R _ IPil 2 |P2| 2 |p 3 | 2 „ mM mjM mmj 

2m + 2M + 2m J U |n-r 2 | U |r 3 -r 2 | \n - r 3 \' [Z> 

with ri, r 2 and r 3 the positions of the FFP, the star and the BP with respect to 
the inertial reference frame x, y and pi, p 2 and p 3 the corresponding momenta. 
The initial position and momentum of the three bodies are: x = 0, y = 0, p x = 0, 
p y = for the star, x = r$ cos </>, y = tq sin <j), p x = —mvo sin </>, p y = mvo cos <fr 
for the BP and x = — 40ro, y = —d, p x = mjv, p y = for the FFP. 

Throughout the numerical experiments we fix the mass ratio of the star - 
planet system equal to the mass ratio of the Sun - Jupiter system and we set the 
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radius of the planet's orbit (r ) equal to the semi-major axis of Jupiter. We use a 
system of units where G = 1, M = 1, ro = 1. It follows that the mass, circular 
velocity and orbital period of the BP are mj — 0.000964583M, vq = 1 and 
T = 2ir respectively The initial velocity of the FFP is determined through cq. 
(TTJ) as a function of d and m and is approximately equal to v = 0.2236«o- The 
equations of motion are integrated through a Roscnbrock scheme, which 
is optimized for stiff ordinary differential equations (Roscnbrock, 1963; Gear, 
1971). 

2.3 Final state — exit tests 

Due to energy exchange between the three bodies the system may end up in 
four qualitatively different states, depending on the total energy of the system 
at the center of mass. In particular there are three types of unbound motion. 
These are: flyby, where the FFP is deflected from the planetary system leaving 
it bound, exchange, where the FFP displaces the BP and takes its place orbiting 
the star and disruption, where the planetary system is destroyed and the three 
bodies move independently. Disruption is possible only when the total energy of 
the system is zero or positive. Bound motion between the three bodies is called 
temporary capture and can occur only when the total energy of the system is 
negative. At temporary capture, the FFP and the BP are temporarily bound 
to the star. Again, due to energy exchange between the bodies at every close 
approach, one of the bodies (BP or FFP) can gain enough energy and escape. 
So, temporary capture can result in indirect flyby (escape of the FFP) or indi- 
rect exchange (escape of the BP). But it is not known a priori if and when a 
temporary capture will "break up" and end in indirect flyby or exchange. 

The state of the system at a given time after scattering is a function of 
the initial configuration, which is determined by the parameters <f> and d, since 
we keep all other parameters fixed. We have to note that the center-of-mass 
energy of the system is a weak function of <j> and d, but it is always negative 
allowing for bound motion (temporary capture) and ruling out disruption. Our 
purpose is to map the initial conditions (</>, d) to the final states of the system 
(flyby, exchange or temporary capture). In order to do this, we choose 500 
values of 4> and 700 values of d uniformly distributed in the intervals < < 1 
(<p in cycles) and — 7?'o < d < 7ro respectively |d ma :r| = 7ro is the maximum 
impact parameter for which a direct exchange can occur and may be estimated, 
as a function of the FFP's mass, using known analytical approximations (e.g. 
see Donnison, 1984a). The equations of motion are integrated for each initial 
conditions set [4>,d) until a certain time ttest- At that time a test is made to 
decide the final state of the system. However, the exit test cannot distinguish 
between direct and indirect, flyby or exchange, since the three final states are 
identified only on the basis of the relative energy of the bodies in pairs at time 
Uest- If En, E23, E13 are the center-of-mass energies of the pairs FFP-star, 
BP-star and FFP-BP, they correspond to the final states as follows: 

• a. Flyby: E X2 > 0,E 2 3 < 0,Ei 3 > 
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• b. Exchange: E 12 < 0, £23 > 0, E 13 > 

• c. Temporary capture: E12 < 0,i?23 < 0, E13 > 

The case of capture of the FFP as a satellite of the BP (E13 < 0), al- 
though possible, never appeared in our numerical experiments. The exit test is 
performed at time t tes t = 80r /0. 2236^0 « 360t*o/«o ~ 57. 3T, where T is the 
unperturbed orbiting period of the BP. At this time the FFP, moving at its 
initial velocity, would have traveled a distance twice its initial distance from the 
star. Typical cases of FFP and BP escapes (since we cannot be sure whether 
temporary capture is a real final outcome, see Section [5]) are given in Figs. [3] 
and [3J It should be noted that temporary capture appears for both prograde and 
retrograde orbits (Gayon, 2009), which, as a rule, originate from corresponding 
initial positive or negative impact parameters. However the opposite is not rare 
(initial negative (positive) impact parameter leading to prograde (retrograde) 
orbits). From this we see that the "hidden" physics below a final outcome (i.e. 
the exact path leading to a certain outcome) is not unique. 

Finally it should be noted that an initial velocity of the FFP 10% higher 
than the "parabolic" one used in the paper (for which the total energy of the 
system remains still negative) results in a decrease of temporary captures. 

3 Scattering is fractal 

3.1 Fractal properties of the initial-value space 

In Fig. [5] we present the final states corresponding to each initial condition in 
(</>, d) for the scattering of a FFP with mass equal to that of the BP (m = mj). 
Each initial condition pair (</>, d) is coded according to the final state of the 
system at the time of the exit test, ttest = 57. 3T. Gray represents flyby, white 
represents temporary capture and black exchange. As it is illustrated in Fig. 
flyby and temporary capture of the FFP are the most probable (> 99%) 
outcomes, while the probability of exchange is much smaller (rj 0.1%). At this 
resolution, it seems that there is a clear boundary between the three basins of 
final states, with the basin of exchange constituting the border between flyby 
and temporary capture. 

Zooming twice in the region defined by the small black box in Fig. [5j we 
see that the boundary between the three basins is not smooth (Fig. [5] bottom). 
Within the region leading to temporary capture (white), there are scattered 
initial conditions that lead to the other two final states, a fact repeated at 
smaller and smaller scales. This indicates uncertainty in predicting the final 
state of the system for regions close to the boundary and a possible fractal 
structure. On the other hand there are regions where the boundary between 
two basins appears to be smooth i.e. the boundary between flyby (gray) and 
temporary capture (white) basins in Fig. [5] bottom-right. The coexistence of 
fractal and smooth (non fractal) boundaries intertwined in arbitrarily fine scales 
is common in Hamiltonian systems (e.g. see Blcher et al., 1990). From the two 
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Figure 3: (Top) An example of FFP escaping case following a period of tem- 
porary capture (indirect - flyby) (FFP: light gray, BP: dark gray, star: black). 
Corresponding evolution of stellarcentric semi-major axis (middle-row) and ec- 
centricity (bottom-row) of the FFP and BP until the FFP's escape. We see that 
the planets jump from one orbital resonance to another (see Fig. 4 as well) 
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Figure 4: (Top) An example of BP escaping case following a period of temporary 
capture of the FFP (indirect - exchange) (FFP: light gray, BP: dark gray, star: 
black). Corresponding evolution of stellarcentric semi-major axis (middle-row) 
and eccentricity (bottom-row) of the FFP and BP until the BP's escape. We 
see that the escape of the BP is the final outcome of a sequence of resonances, 
which drive the FFP to lower values of (a, e), while the BP's elements increase. 
Note: The upper left corner panel does not depict the escape of the BP, which 
happens at a much later time. 
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Figure 5: (Top) Initial- value space for the scattering of a FFP with mass 
m = mj. Each initial condition cj>,d is coded to represent the state of the 
system after the scattering, at time ttest = 57. 3T. By this time, the majority of 
initial conditions have led to either flyby (gray) or temporary capture (white). 
Initial conditions corresponding to exchange (black) constitute the boundary 
between the other two basins. Note that temporary capture (white) is not a 
definite outcome and the corresponding orbits may later on decay to either flyby 
or exchange. (Bottom-left) Magnification of the region defined by the black box 
and (bottom-right) a further magnification of the sub-region contained in the 
white box. The occurrence of scattered initial condition leading to flyby and 
exchange within the basin of temporary capture (white), leads to great uncer- 
tainty in predicting the final state of the system for initial conditions close to 
the boundary. The two successive magnifications show an apparent self simi- 
lar structure of the boundary, with fractal boundaries intertwined with smooth 
boundaries at finer and finer scales. 
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successive enlargements shown in Fig. [5] it is also apparent that the non-smooth 
border looks qualitatively self-similar for at least two scales of magnification. 

3.2 Estimating the uncertainty exponent 

In order to estimate the boundary (between basins) dimension, we apply the 
final state sensitivity method introduced by Blcher et al. (1990). The basic idea 
behind this method is to estimate the probability, p(s), to make a wrong pre- 
diction of the final state when there is an uncertainty e in the initial conditions. 
If the boundary is a smooth 1-D geometrical object (i.e. a line) this probability 
would be exactly proportional to e. In a different case (i.e. a fractal boundary) 
we expect this probability to scale differently with e. In any case, the scaling 
of p with e gives the co-dimension of the boundary, a, from which the capacity 
dimension of the boundary follows as d = D — a, with D the dimension of the 
embedding space (see Bleher et al. 1990). The method is the following: 

In a box [<fi, d] on the initial-value space we select at random 25,000 "central" 
initial conditions (0, d). Each "central" initial condition is perturbed along the 
d axis by an amount e, so that we get two "side" initial conditions (</>, d + e) and 
(<fi, d — e). If the three initial conditions (the "central" and the two "side" ones) 
do not all lead to the same final state, the "central" condition is considered as 
e-uncertain. For every e we calculate the fraction, /(e), of e-uncertain "central" 
initial conditions and plot it as a function of e in log-log axes. The points 
(e, /(e) ) should lie on a straight line log/ = aloge + 6, where the slope, a, 
is the co-dimension of the boundary, named in this case uncertainty exponent. 
The capacity dimension, d, of the boundary is d = D — a, with D = 2 (the 
dimension of the initial- value space). A smooth boundary has a = 1, so that 
d = 1 (a line). If a < 1 (as we find in most cases here), the boundary is fractal 
and d > 1. 

To estimate the uncertainty dimension we select 10 sub-regions of the initial- 
value space, shown in Fig. [SJ In each region we compute the capacity dimension 
of the boundaries of the basins of (i) exchange and (ii) flyby, considering in each 
case the remaining two basins as a single one. The scaling of / with e for region 
3 of Fig. |6] (left) and for the basin of flyby is shown in Fig. [6] (right). For 
e < 10 -5 the function /(e) is very well fitted by a straight line for at least 3 
orders of magnitude of e, which is a strong indication that our results show a 
real fractality of the boundary and are not an artifact of the method. The fit 
gives a = 0.257 ± 0.014 and, therefore, d ~ 1.743. The corresponding values 
of the uncertainty exponent in the other regions of Fig. 6 range in the interval 
0.209 < a < 0.346. 

It should be noted that the uncertainty exponent, a, is a decreasing function 
of the test time, t test . This is because, as we will see in more detail below, 
temporary captures decay over time in an exponential fashion. Thus, on much 
longer time scales, there appear to be only two possible outcomes - instead of 
three - and the boundary becomes smoother. However, as we will see in the 
following, the number of temporary captures does not seem to go to zero on 
reasonably long time scales (see section [5]). 
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Figure 6: (Left) Black boxes denote 10 regions of the initial value-space for 
which the uncertainty exponent was estimated (dimensions of the regions 8(j) 
= 0.003 and Sd = 0.06). (Right) Fraction of uncertain initial conditions, /. 
as function of the uncertainty, e, for region 3 and for the boundary of flyby, 
considering the other two basins as one. The uncertainty exponent is estimated 
to a = 0.257 ± 0.014, which gives a boundary dimension of d « 1.743. 

4 Final state statistics 

4.1 Probabilities of exchange, flyby and temporary cap- 
ture 

We have performed a large number of numerical experiments for various values 
of the ratio m/mj and in this section we give the statistics at time t tes t = 57. 3T, 
concerning (direct) exchange, (direct) flyby and temporary capture as functions 
of m/mj. 

The probabilities of exchange and flyby as functions of the FFP's mass are 
shown in Fig. [7] Each point in the graphs represent a simulation of 350,000 
initial conditions. Roughly half of the cases end up in direct flyby, while a small 
percentage (of the order of 1%) ends up in direct exchange. The rest of the 
cases (s=s 50%) correspond to temporary capture of the FFP. It is interesting to 
note that, for a relatively large interval, the exchange and flyby probabilities 
depend linearly on the FFP's mass. Another interesting point is that there is a 
lower mass limit (m — 0.5m j), below which direct exchange cannot occur. An 
explanation of this fact is given in section (|4.3[) . A similar result concerning the 
lower mass limit for exchange has been reported by Boyd and MacMillan (1993) 
in the case of scattering of a star from a binary star of equal mass members, 
within a configuration similar to the one used here. In particular, these authors 
found that the probability of exchange drops to zero for masses of the incoming 
body between 0.3 and 0.8 m s t a r (M). For the problem of equal binary masses 
the probability of exchange increases rapidly with the mass of the incoming 
body, since in this case either one of the binary members can escape. In our 
case, the masses of the BP and the FFP are much smaller than that of the 
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Figure 7: Fraction of orbits that have resulted in flyby and exchange by the 
time of the exit test (t tes t = 57. 3T) as functions of the FFP's mass. The third 
outcome, namely temporary capture of the FFP, is the supplement of the other 
two sets and is an increasing function of the FFP's mass in this range of m. 
Note that the probability of exchange falls to zero for m < 0.5mj (see section 
1431) . 



star and the probabilities of exchange, flyby and temporary capture do not vary 
strongly with the mass of the FFP. 

4.2 Orbital distribution in temporary capture 

In temporary captures the orbital elements of both BP and FFP vary within a 
wide range. Shortly after scattering, the majority of captured FFPs have very 
elongated orbits (e > 0.99) with semi-major axis 1.998ro < a < 3.3 ■ 10 7 ro, 
while the BP is slightly perturbed form its initial orbit, a result corroborated 
by Perets and Kouwcnhoven (2012). In 90% of the cases the orbits of the two 
planets do intersect. However, there are also a few cases where both planets 
have moderate values of a and e. Fig. [5] shows such cases for two values of 
the FFP's mass. Depending on the mass, these orbits make up from 0.01% 
(??i = 2mj) to 0.2% (m = 0.5m,/) of the initial sample of orbits. Finally, 
although the eccentricity and semi-major axis show important variations, there 
is an approximate conservation of pericenter distance for both bodies, with 
q « 1.0075 for the FFP and q « 1.1057 for the BP. 

As time passes (t = 600, 000T) the picture changes, since now 75% of the 
orbits do not intersect, as most of the interacting orbits have decayed to flyby 
or exchange. In all these orbits the FFP has pericenter distance in the range 
2?'o < q < 1291ro, while the BP has semi-major axis in the range 0.9ro < a < 
1.2ro and moderate eccentricity. There are also a few orbits where the BP and 
FFP exchange positions (the FFP jumps to an inner orbit and the BP goes 
outwards). As we show in section [5l there exist bound orbits that survive for 
very long times. 

It should be noted that tidal interactions between the FFP and the star, in 
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Figure 8: Stellarcentric orbital elements of the FFP (black) and the BP (gray) 
in temporary capture, for orbits for which both planets have moderate values 
of semi-major axis and eccentricity (a < 50ro, e < 0.7). The elements are 
computed at time ttest = 57. 3T. 
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Figure 9: Eccentricity (e) as a function of the semi-major axis (a) of the FFP 
in the event of direct exchange, for four different values of the FFP's mass: 
m = 0.5toj, m = lmj, m = 2mj and m = 8mj. The stellarcentric osculating 
values of a and e are computed at time ttest = 57. 3T for the orbits that have 
resulted in exchange by that time. The solid lines are constant pericenter and 
apocenter curves e = 1 — 1/a for a > 1 and e = l/o — 1 for a < 1). 
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the case of extended bodies, are expected to stabilize the FFP's orbit, if the 
percenter distance, g, is relatively short. Due to this process, named tidal circu- 
larization, the eccentricity of a planet is decreased through the combined effects 
of energy dissipation, through tidal friction, and conservation of angular mo- 
mentum. From analysis of our results we have found that, after t = 200, 000T, 
6% of the orbits have q < 0.1 AU; this fraction drops to 2% after t = 600, 000T. 
The fraction of orbits with small pericenter distance is of comparable order of 
magnitude to the long-time surviving orbits found in section [5] and, therefore, 
are not negligible in the final analysis, since they might increase the fraction of 
FFPs eventually captured. 

4.3 e vs. a in direct exchange 

Fig. [5] shows an interesting "regularity" in the event of exchange. The eccen- 
tricity of the captured FFP is plotted as a function of its semi-major axis, for 
various values of the FFP's mass. In all cases a and e are contained in a region 
in a, e space bounded by two curves, e = 1 — (1/a) for a > 1 and e = (1/a) — 1 
for a < 1. The first one corresponds to constant pericenter distance q = 1.0 
and the second one to constant apocenter distance, a(l + e) = 1. By setting 
e = 1 in the second equation we obtain a lower bound for the semi-major axis 
of the FFP, which equals to a = 0.5ro, independent of the FFP's mass. There is 
also a less clear upper bound of the value of the semi-major axis, which seems 
to depend on the FFP's mass. This upper bound is a result of energy conser- 
vation. In particular, the energy taken by the escaping BP is gained at the 
expense of the star-FFP energy. Since the binding energy of the initial star-BP 
pair is Ej = —(GMmj)/2ro and the FFP has initially zero binding energy, 
the binding energy of the new pair will be equal to or less than Ej, so that 
-(GMm)/2a < -(GMmj)/2r a . It follows that the semi-major axis of the FFP 
will be a < (m/mj)ro- According to this, an FFP with m = mj will have at 
most a = lro, which is apparently the case in Fig. [5] Now consider an FFP 
of mass m < 0.5mj. This will have semi- major axis a < 0.5ro. However this 
value is below the lower bound a = 0.5tq (imposed by the lower bounding curve 
(e = (1/a) — 1)). Therefore a direct exchange cannot occur for masses of the 
FFP below 0.5m j. 

An interesting point, coming out from Fig. 9, is that in some direct exchange 
outcomes for m/mj = 2 and m/mj = 8 we observe that a > (m/mj)ro. We 
have looked in detail onto this and we have found that the additional energy 
AE of the FFP-star pair is taken from the FFP-BP one, which is non-negligible 
when the FFP is relatively massive. In Fig. [10] we see that AE obeys an 
approximate relation of the form AE oc 1/m. 

5 Exponential decay of temporary captures 

A temporary capture configuration at t tes t can decay through the escape of 
either the BP or the FFP. In order to study the decay process, we continued the 
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Figure 10: Energy loss (per unit mass) from the FFP-star pair as a function 
of the FFP's mass. We see that a curve of the form AE « 1/m is fitting the 
data with a coefficient of determination R 2 ~ 0.9. AE is computed for the 
orbits with semi-major axis that exceeds the upper bound a < m/mj)ro- AE 
is the difference between the true binding energy of the FFP and the energy 
that corresponds to the upper bound. 
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Figure 11: (Left) Fraction of orbits that decay to exchange vs. time. (Right) 
Fraction of orbits that decay to flyby vs. time. N = 6, 500 is the initial sample of 
orbits that were in temporary capture at time t tes t = 57. 3T and were integrated 
up to time t = 600, 000T. 



numerical integration of 6,500 cases that were in temporary capture at ttest = 
57. 3T up to t final = 600,000T. Fig. [TT] shows the fraction of orbits that have 
resulted in indirect exchange (left) and flyby (right) as functions of time. We 
see that more than 90% of the orbits, that were in temporary capture at time 
ttest j will result in flyby and approximately 4% of them will result in exchange 
as t —> oo. Fig. [T2l shows a lin-log plot of the fraction N(t)/N of orbits, that are 
still in temporary capture after time t, vs. t. We see that, after t = 200, 000T, 
there seems to be an exponential fall-off of this fraction. By fitting a least square 
curve of the type N(t)/N — No/N = (3 ■ e~ xt to the data we find (through a 
X 2 test, see Fig. Q2) that the best fit is: N(t)/N - 0.016 = (3 ■ e" 410 " 6 *. This 
means that there is a fraction of orbits, N /N ~ 0.016, that are expected to 
remain in temporary capture for very long times. It should be noted that the 
fractal properties of the basin boundaries persist, even after the decay of the 
temporary capture configurations to either indirect flybys or indirect exchanges. 



6 Prograde and retrograde orbits 

Fig. Q2] is a plot of the initial (0, d) values that correspond to prograde (gray) 
and retrograde (black) orbits of the two planets, that remain bound up to 
t = 60, 000T. As expected, orbits with positive impact parameter end up mostly 
as prograde and those with negative impact parameter mostly as retrograde. 
However, there are orbits for which the FFP or the BP change their initial revo- 
lution direction. This happens for low values of the impact parameter, between 
— 2tq and 2ro for which the FFP's orbit passes very close by the BP's orbit (see 
Fig. Q2] where we have prograde orbits (gray) with d < and retrograde orbits 
(black) with d > 0). In the majority of the cases it is the FFP that changes 
its revolution direction. This is to be expected, since the FFP orbit is nearly 
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Figure 12: (Left) Fraction of orbits that remain in temporary capture until 
t = 600, 000T and the curve that fits best the data. (Right) Coefficient of 
determination, R 2 , as a function of Nq/N. The best fit occurs for Nq/N = 0.016, 
where Nq/N is the fraction of orbits that are expected to remain in capture as 
t — > oo. N = 6, 500 is the initial sample of orbits that were in temporary capture 
at time t tes t = 57. 3T and were integrated up to time t = 600, 000T. 
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Figure 13: Initial conditions corresponding to orbits that remain in temporary 
capture until t = 60, 000 T. Grey indicates prograde orbits and black indicates 
retrograde ones. 
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parabolic and therefore can easily become a very eccentric ellipse. The angular 
momentum exchange between the FFP and the BP depends on the minimal 
distance and relative phase of the two bodies, during their encounter. Thus, the 
smaller the impact parameter, the closer the two bodies can get. This increases 
the magnitude of the angular momentum exchanged and - depending on the 
relative phase - increases the probability that the FFP's orbit will change its 
revolution direction around the host star. 

7 Summary and conclusions 

In the present paper we have studied the gravitational scattering of a FFP by 
a star-planet pair and we have shown that it is chaotic, in the sense that the 
boundary of a region in initial conditions that leads to a specific outcome (i.e. 
flyby, capture or exchange) is fractal. In this way our statistical approach is 
fully justified. 

One interesting result, coming out from our analysis is that the probability of 
exchange, flyby or temporary capture is, for values of the FFP's mass spanning 
one order of magnitude, a linear function of the FFP's mass, decreasing in the 
case of flyby and increasing in the cases of exchange and temporary capture. 
It should be noted that for masses of the FFP below 0.5mj exchange is not 
possible, due to energy conservation. 

Another interesting result is that, in the case of exchange or temporary cap- 
ture, FFPs can be captured in orbits with moderate values of (a, e), with the 
simultaneous excitation of the BP's orbital elements. In the case of temporary 
capture, orbits where both FFP and BP have moderate values of (a, e) make 
up w 0.2% of all captured orbits (depending on the FFP's mass). Restricting 
ourselves to non-Hill stable interactions, there is a cumulative 4% probability 
that a FFP, passing within a distance of 7tq from the star (where ro is the radius 
of the BP's circular orbit) will be captured in orbit around the star. This cumu- 
lative probability is the sum of direct exchange (0.2%), indirect exchange (2%) 
and asymptotic captures (1.6%). According to Donnison (2008), the extension 
to 3D should increase the probability of capture of the FFP by the star-BP pair. 

The scattering induces an excitation of the orbital elements of the final 
system in both exchange and temporary capture. In exchange, the FFP can 
be captured in highly eccentric orbits (e > 0.9) with semi-major axis up to 
a = (m/mj)r (where m the mass of the FFP and mj, ro the mass and semi- 
major axis of the bound planet). In temporary capture, the majority of the 
FFP orbits that survive up to t « 600, 000 orbital periods of the BP have very 
large and elongated orbits (e > 0.97 and a > 5.000r ) with the BP being only 
slightly perturbed from its initial orbit. Nevertheless 75% of these orbits do not 
intersect. 

The above results come out from a model with two main simplifications: (i) 
co-planar configuration and (ii) absence of dissipation and some minor simplifi- 
cations: (iii) circular orbit of the BP, (iv) fixed planet to star mass ratio and (v) 
parabolic orbit of the FFP. It would be of great interest to perform further "ex- 
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periments" , by introducing (i) the third dimension and (ii) a mechanism of tidal 
dissipation acting after the capture, in order to study the long-time evolution 
of FFP's in high-e and high-a orbits. 
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